Single file diffusion in macroscopic Wigner rings 
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The single file diffusion in a circular channel of millimetric charged balls is studied. The evolution 
in time of the mean square displacement is shown to be subdifTusive, but slower than the power-like 
t^^'^ behavior observed in circular colloidal systems or predicted in one-dimensional infinite systems. 
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I. INTRODUCTION 

Diffusion of particles in a single channel where the mu- 
tual exchanges are forbidden is known as single file dif- 
fusion (SFD). This physical situation is encountered in 
various fields, for instance ID hopping conductivity, ion 
transport in biological membranes or molecules channel- 
ing in zeolithes. In these geometries with low dimension- 
ality, the particles in repulsive interaction are not able 
to cross each other, and the particles remain correlated 
even at long times. 

Such systems were generally modeled by an infinite 
set of particles with hard interaction and diffusing on a 
line 0, Q : the mean square displacement Ax^ has been 
proved to grow at long times as t^/'^. More recently, it 
was shown that this behavior is also encountered in over- 
damped systems with arbitrary repulsive interactions, 
providing the correlation length between the particles is 
of finite range Q- On the other hand, although the main 
theoretical results for single file diffusion were obtained 
many years ago, experiments displaying such subdiffusivc 
behavior are lacking and the obtained results are often 
conflicting. For instance, NMR studies j4u5i] and quasi 
elastic neutron scattering experiments [J, Q on organic 
molecules in porous materials either conclude to subdif- 
fusivc transport 0, Q or to classical diffusion ^ JJ , with 
apparently the same experimental conditions. More re- 
cently, the movement of colloidal particles in a circular 
channel has been studied. This channel is obtained by 
photolithography 0, Q or optical tweezer circular trap- 
ping Q. While classical diffusion is always observed for 
short times, the long time behavior is not as well clearly 
analyzed. Some authors suggest it grows as 

^1/2 like the- 
oretical SFD processes for infinite systems 0, 0] whereas 
others explain this subdiffusion by hydrodynamic cou- 
pling effects 8]. On the other hand, the theoretical study 
made in Ref. 2] for a ID line of hard-core interacting par- 
ticles with periodic boundary conditions shows that the 
mean square displacement grows linearly, as for a free 
diffusion. 

In order to rule out these ambiguities, we present here 
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FIG. 1: Photography of the experimental circular channel 
with A'' = 12 balls inside. 



the diffusion of macroscopic charged metallic balls elec- 
trostatically interacting and moving in a circular chan- 
nel whose width forbids any crossing, while a mechani- 
cal shaking induces an effective temperature. This sys- 
tem presents many advantages : it gives the opportu- 
nity to study the diffusion at very long times, it allows 
to suppress any hydrodynamic effect and, if necessary, 
to tune the interacting forces. In this experiment, the 
inter-particle interaction is similar to inter-vortices 
interaction in superconductors jlQ] . This interaction dif- 
fers from the dipolar interaction of colloids or from the 
hard interaction introduced in many theories, thus we can 
evaluate the influence of the interaction characteristics on 
the SFD behaviors. Note that no hard-core collisions are 
observed in our system. 

Our main result is to exhibit subdiffusive behavior 
slower than the t^/^ behavior predicted by theory and 
mentioned in colloidal systems. 

In Sec.niwe present the experimental set-up and prove 
its ability to describe diffusion processes by studying two 
basic cases. To identify the effects of the circular confine- 
ment on the diffusion in order to distinguish them from 
the inter-particle interaction contributions and to deter- 
mine the best conditions to obtain the required ID move- 
ment, we have studied the diffusion of a single ball in the 
channel. The obtained results are reported in Sec. IIIII 
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Sec. IIVI is devoted to the diffusion of = 12 and 16 
interacting particles and to the comparison with the be- 
haviors previously presented in literature, in particular 
those obtained with colloids. 



II. EXPERIMENT VALIDATION 

In this experiment, millimetric stainless steel balls (of 
radius R = 0.4 mm and weight m = 2.15 mg) are located 
on the doped-silicon bottom electrode of an horizontal 
plane condenser (of height h — 1.5 mm). A metallic 
frame intercalated between the two electrodes and in con- 
tact with the bottom one confines the balls in a circular 
channel. Its external diameter and its width are respec- 
tively 10 mm and 2 mm (Fig^). This forbids any crossing 
between the balls. In order to charge the balls, a tunable 
voltage Vc of about IkV is applied to the top electrode. 
We then get a system where the interaction between the 
balls as well as the confining potential are easily ad- 
justed. The interaction potential V{r) between two balls 
has been shown to be well described through a modified 
Bessel function of the second kind : V{r) — AKo{r/X) 
where the screening length A is about 0.3h — 0.5 mm and 
independent from Vc, as discussed in details in Ref. [T(tI |. 
Direct capture of the position of the balls is made by a 
camera placed above the experimental device. The typ- 
ical time between two snapshots is selected between 15 
ms and 150 ms and series of 10,000 images were recorded. 
The individual trajectories of these particles are thus di- 
rectly determined over very long times (comparing with 
the typical relaxation time that will be determined in the 
following) . 

To introduce the thermal noise, the whole cell is 
fixed on loudspeakers supplied by a white noise voltage. 
Thanks to the friction between the bottom electrode and 
the balls, a erratic and non spatially correlated move- 
ment is conveyed to the latter. It has been previously 
shown that the ensemble of balls obey Boltzmann statis- 
tics, where the shaking amplitude stands for the temper- 
ature 11]. The effective temperature is thus determined 
in situ by measuring the mean square position of a single 
ball confined by a circular frame. 

The diffusion of the particles is characterized by the 
evolution with time of their mean square displacements 
(m.s.d.) given by, if the considered coordinate is x : 



(1) 



where the brackets ( ) denote ensemble averaging. In 
practice, we get from the trajectory of one ball a set of 
statistically independent trajectories by shifting the time 
origin. We emphasize that such a method implies un- 
avoidable averaging over the distribution of initial con- 
ditions, which can lead to mean square displacements 
slightly different than the usual ones (see the appendix). 
Furthermore, when considering the orthoradial coordi- 
nate 6 (in radians) we will consider the cumulated an- 
gle and not the modulo 2t: angle, since we want to con- 
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FIG. 2: (a) Typical mean square displacement of a free ball 
(log scale) ; (b) typical mean square displacement of a ball 
trapped in a well constituted by a circular frame of radius 5 
mm. The applied voltage Vc is 1000 V. 



sider the ring as an easy realization of a ID line. As in 
Ref. [1, ^, , the total orthoradial displacement of a ball 
is taken into account, without any substraction of a mean 
displacement of the system. 

In order to evaluate the ability of our experimental set- 
up to explore diffusion processes, two simple cases for 
which the diffusion behaviors are well established have 
been studied. 

The free diffusion of a ball moving on the bottom elec- 
trode without any confinement and any applied electro- 
static potential was first studied. Whatever the effec- 
tive temperature, the m.s.d. Ax'^(t), x being the ball 
abscissa, presents the same typical time variation, pre- 
sented on Fig. a). This m.s.d. varies linearly with time 
at long times, whereas the short time behavior exhibits 
a power increase with 1 < /3 < 2. These behaviors are 
in agreement with a free diffusion described by the clas- 
sical Langevin equation which predicts Aa;^(t) ~ 2Dt 
at long times where D is the diffusion constant. Since 
the temperature is independently measured, the effec- 
tive damping coefficient 7 can be determined from the 
fiuctuation-dissipation theorem m'^D = fc^T. The re- 
lated relaxation time — I/7 which marks the crossover 
between the short and long time regimes is found to be 
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equal to about 100 ms, which is coherent with a rough 
estimation on the curve of Fig.|21[a). 

This value which is quite close to the smallest snapshot 
time explains why the theoretical behavior in predicted 
for very short times {t tr = 1/ 7) may not be observed, 
the behavior observed at short times corresponding actu- 
ally to the progressive transition between a regime to 
a linear one (this analyze is confirmed by the measured 
mean square velocity which is always constant in time as 
expected in the Langevin theory for t > tr). 

The second test was to follow the diffusion of a particle 
confined in a well created by a circular frame (of radius 
5 mm). In such a case, the ball moves in a parabolic 
confining potential whose stiffness K can be determined 
from the radial distribution of the ball Fig. EKb) 

presents a typical radial m.s.d. of a ball in such potential. 
After a rapid increase, it reaches a constant value which 
varies linearly with the effective temperature as expected. 
These behaviors are those predicted in a Langevin for- 
malism in which a parabolic potential models the con- 
finement. At short times {t < (if/m)^/^), the ball begins 
to follow a free diffusion since the particle has no time 
to explore the whole confining potential whereas at long 
times (t > (if/m)^/^) the experimental m.s.d. is depen- 
dent on the confinement and reaches 2kBT/K (see the 
appendix). In the same time, the angular m.s.d grows 
linearly with time, which characterizes a free diffusion 
along the orthoradial direction. 

The accordance between the theoretical behavior and 
the observed one confirms the ability of our experimen- 
tal set-up to give relevant and self coherent parameters 
about particle diffusion and that the Langevin formalism 
describing Brownian motion can be used. 



III. DIFFUSION OF A SINGLE BALL IN A 
CIRCULAR CHANNEL. 

In our experiments, the interacting balls move in a 
electrostatic confining potential which looks like a circu- 
lar gutter. In order to distinguish possible effects due 
to this channel geometry from those resulting from the 
inter-particle interaction, we have preliminary studied 
the diffusion of a single ball in such a circular channel. In 
particular, the radial confining potential has been stud- 
ied for different applied electrostatic potentials in order 
to get the best opportunities for a quasi ID movement. 
In the same time, the coupling of the small radial dis- 
placement with the orthoradial movement inherent to the 
circular movement has also been evaluated. 

At low temperatures, the ball oscillates radially with 
a slow orthoradial drift. As the temperature increases, 
the orthoradial displacements grow, the ball runs over 
the whole channel with long free jumps between two se- 
quences of radial oscillations. Note that no bounces on 
the channel walls have been observed whatever the tem- 
perature, the ball is only confined by the electrostatic 
potential. Its typical shape, determined from the distri- 
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FIG. 3: Confining radial potential in the circular channel ob- 
tained from the radial distribution of the balls, for two dif- 
ferent numbers of balls and two different applied voltage Vc 
: 1 ball, T4 = 800 V (•); 1 ball, K = 1000 V (a); 12 balls, 
Vc = 800 V (o); 12 balls, K = 1000 V (a). r is measured 
relatively to the mean radial position. 



<l 0.04 



1 ' 1 ' 1 


i 1 i 1 i 1 1 1 1 1 i 1 i 




i^*^ ^""^ V=600V 

C 








^ V =800 V 

c 




V=1000 V " 

c 


I.I.I 


,1,1.1. 



200 400 600 800 1000 1200 1400 1600 1800 



t (ms) 

FIG. 4: Radial m.s.d. for a single ball in the circular channel 
for three different applied voltages. Solid lines show the fits 
to equation 121 



bution of radial positions, is presented in Fig. What- 
ever Vc, it is parabolic with a stiffness K that increases 
with the applied voltage while the mean radial position 
remains unchanged. Let us indicate that for too small Vc, 
the electrostatic confinement looks more like the flat geo- 
metrical channel profile than a parabolic well. Therefore, 
we will work in the following with high Vc (— lOOOV) in 
order to reinforce the ID movement. 

For a ID parabolic potential, Langevin equation can 
be analytically solved and its characteristic coefficients 
determined by fits with the experimental data. The evo- 
lution with time of the radial m.s.d. of a single ball in the 
channel is shown on Fig. ^ At short times, it presents 
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a rapid increase followed by small damped oscillations 
and finally reaches a constant value at long times. It 
is very well fitted by the theoretical expression obtained 
through Langevin formalism in the case of low friction 
(see the appendix), 

Ar^it) = 2{{rl)) [l - e-^*/2( cos(^t) + ^ sin^))] , (2) 

where ((rg)) = ksT/K is the mean square starting 
position, 7 the damping coefficient and uj = a/wq — 7^/4 
(with Wo = K/m) the effective frequency of the oscilla- 
tions of the ball. 

As expected, the larger the applied potential Vc is, 
the better the fits are. The frequency w, which is 
temperature-independent and equal to about 20 Hz 
within the explored Vc range, is in accordance with the 
independently measured frequency obtained from the 
Fourier transform of the ball trajectory, which confirms 
the validity of the use of the Langevin equation. The K 
values determined from these frequencies are in perfect 
agreement with those evaluated from the radial distri- 
bution variance and varies as as expected from the 
determination of the interaction between a ball and the 
wall [13. 

From the fits, the damping coefficient 7 is found to be 
independent from Vc and equal to about 10 s~^. This 
value is of the same order as the one for free diffusion, 
which indicates that the relaxation time r^j = 7"^ results 
uniquely in the case of a unique ball from the friction 
process on the bottom electrode. Within the temperature 
range, 7 is slightly temperature-dependent but, because 
of the respective range of and 7^/4, this dependence 
is not significant for w. 

Note that the main feature of the presented fits is that 
the time t/j, which is an important time of reference when 
analyzing the orthoradial movement, can be precisely de- 
termined for each experiment. 

Another point is the coupling between the orthoradial 
and the radial movement in the circular channel. The 
variations of the angular m.s.d. for Vc — 1000 V is pre- 
sented in Fig. 13 The long time [t ^ tb) behavior is 
linear, which indicates a free orthoradial diffusion. This 
behavior is observed whatever the temperature. Note 
the oscillations of the radial m.s.d are not observed here, 
which strongly evidences that the orthoradial and radial 
movements are independent. This was suggested by the 
independence of K with the effective temperature, be- 
cause otherwise the effective stiffness would depend on 
the temperature through the orthoradial velocity. More- 
over this angular free diffusion time dependence unam- 
biguously confirms that the ball is only confined by the 
electrostatic potential and that the balls are not bounc- 
ing on the walls, otherwise different dependencies would 
be found 

Finally, from the measurement of the diffusion con- 
stant of the orthoradial movement and the measure- 
ment of the damping coefficient through the independent 
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FIG. 5: Radial m.s.d of a ball in a 12-balI ring, with Vc = 1000 
V and T = 3.5 x lO" K (■), 7.5 x 10" K (□), 11.5 x lO" K 
(•), 17.5 X 10^^ K (o). Solid lines shows the fit to equation|5| 

study of the radial movement, one can check that the 
fluctuation-dissipation relation is obeyed, which proves 
that the thermal equilibrium is reached (as also shown 
in Ref. jllj) and justifies the use of the latter relation to 
determine any missing parameters in a given experiment. 

Thus the movement of one particle trapped in the cir- 
cular channel is twofold. Along the radial direction it is a 
ID diffusion in a parabolic trap whereas it is a free angu- 
lar diffusion along the channel. Moreover, this indicates 
that the formal coupling which exists between these two 
coordinates in the Langevin equation, is not relevant in 
our experiment. So we can use this set-up in order to 
study SFD processes through the orthoradial diffusion of 
N interacting particles. 

IV. DIFFUSION OF INTERACTING 
PARTICLES IN A CIRCULAR CHANNEL 

We now consider rings oi N = 12 and TV = 16 interact- 
ing particles. These two systems were selected because 
in our channel geometry, these rings remain circular and 
do not present "zig-zag" configurations which could al- 
ter the SFD effects [T3 . Furthermore, these rings also 
appear as outer shells in Wigner islands whose dynamics 
will be presented in a next paper pl| . 

The typical radial confining potentials for iV = 12 are 
shown on Fig. |3| As for = 1, they are parabolic and 
temperature independent and their stiffness increases as 
Vc - However we observe a stiffening with N, which is in- 
dependently confirmed by the increase of the frequencies 
obtained from the fit of the radial m.s.d. by the expres- 
sion |21 (Fig. [SJ. The frequencies ~ o^o are about 30 
s^^ for A^ = 12 and 16 whereas they are about 20 s"""^ 
for A^ = 1. This increase of the frequencies with N for a 
small number of balls is in accordance with the numerical 
study presented in Ref. jl3- Consequently, the ID char- 
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T (10 K) 

FIG. 6: Evolution with the temperature of the damping co- 
efficient 7 for Ai" = 12 balls. Solid line shows the fit to the 
Guzman- Andrade law. 




FIG. 7: Angular m.s.d. in log scale for a single ball (□, Vc — 
1000 V, T = 15.1 X 10" K) and a baU embedded in a 12- 
ball ring (solid line, K = 1000 V, T = 3.5 x lO" K). 6 is in 
radians. 



acter of the angular movement is amplified in the case 
of N interacting balls. Moreover, this stiffening makes 
the parabolic modelization of the radial well all the more 
exact, and thus the determination of the damping co- 
efficient 7 through the radial m.s.d. fits all the more 
accurate. As shown on Fig. |^ 7, which remains inde- 
pendent of the applied potential Vc and of the number 
of balls, is a decreasing function of the temperature. Let 
us indicate that this last dependence, well fitted with 
the law 7 = Ae^^^ suggests that our "thermal bath" 
looks more like a liquid than a gas for which 7 increases 
with T We underline it does not imply that the 

electrostatic coupling between the balls is modulated by 
hydrodynamic effects as in Ref. j^]. 

Along the orthoradial direction, the distribution of the 
position of each ball is gaussian with a width propor- 
tional to The amplitude of the displacement is 
small (about 2 mm) , so the curvature of the channel is not 
relevant. The variation with time of the orthoradial m.s.d 
is shown on Fig.[7|for N = 12. As for the single-ball case, 
no oscillations are observed, confirming the complete de- 
coupling between the radial and orthoradial movement. 
After a t'^ increase at short times, a fast transition to 
a subdiffusive behavior at long times is observed. How- 
ever, and it is the strongest difference with the results 
presented for colloidal systems, this increase with time 
is slower than t^/^ at long times. The same behavior is 
observed whatever N and Vc- This continuous diffusion 
slowing down has never been mentioned previously, all 
the experiments on colloids concluding in the existence 
of only two distinct regimes : free diffusion characterized 
by a linear behavior at short times and a t^/^ behavior 
at long times, attributed either to SFD or to hydrody- 
namics coupling effects. We can note also that, contrary 
to the studies on colloids, no clear linear regime is ob- 
served in our experiment. Indeed, the direct interaction 



time which marks the switch between the normal and 
the subdiffusive regimes and qualitatively corresponds to 
the time when a particle "feels" the moving well around 
itself is quite short and of the same order as th which 
characterizes the crossover between the P and the linear 
free diffusion regimes. 

Different reasons could explain these differences. The 
first one results from the characteristic times of the dif- 
fusion itself. In the colloidal systems, the direct inter- 
action time (about 10 s) appears on the reported data 
much larger than (which should be lower than 100 
ms), and then the linear regime can appear and be ob- 
served during one decade as in Ref. 9]. By contrast, in 
our case, the interaction is probably more long range, 
thus the presence of adjacent particles is more rapidly 
efficient and the crossover with the subdiffusive regime 
appears sooner. This assumption is confirmed by the 
decrease of the crossing time as the number of balls in- 
creases, in accordance with the enhancement of the re- 
pulsive interaction. Thus the linear regime is drastically 
reduced, hidden by the rapid emergence of the subdiffu- 
sive behavior. The second reason may be the limitation 
in time in the colloids experiments. In the latter, the 
subdiffusive behavior is put in evidence during another 
decade after the direct interaction time and is fitted as 
t^^'^. However, the m.s.d. time evolution in the figures 
presented Ref. 0, exhibits ultimate points under the 
t^/^ fit curves which suggests a slight slowing down be- 
yond. This slowing down could be observed in our exper- 
iments and not in colloidal systems simply because they 
are performed up to a time of three order of magnitude 
higher than our direct interaction time which is roughly 
equal to Tji ~ 100 ms. Longer recording in colloidal sys- 
tems could give behaviors similar to ours. 

Our results show unambiguously that the behaviors ob- 
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served for colloidal systems are not due to hydrodynamic 
effects but describe real SFD. However, the fundamental 
origin of the observed slowing down is not clear yet. The 
theoretical results which predict t^/"^ at long times focus 
on infinite systems and have been performed either for 
hard-core interactions or for overdamped systems. Nei- 
ther our Wigner rings (the oscillations observed in our 
radial m.s.d. prove moreover that our damping is weak) 
nor the colloidal systems correspond to these assump- 
tions. So the usual theoretical predictions are not fully 
adapted to analyze our results which may not be dis- 
cussed in this simple frame. We suggest that the observed 
slowing down could result from the periodicity of the sys- 
tem associated with an interaction that would be neither 
hard-core nor linear. Indeed, it has been proved that for 
hard-core interacting particles, the diffusion is normal for 
periodic systems On the other hand, for a linear in- 
teraction, the system of coupled Langevin equations can 
be solved in normal coordinates, and the contribution of 
the mode corresponding to the collective rotation gives 
a classical linear diffusion. Since neither our system nor 
the colloidal ones exhibit such a diffusion, this suggests 
that the distance dependence of the interaction is of high 
importance for periodic systems. These situations could 
introduce specific correlations between particles, not de- 
scribed in the usual models. It would be interesting to 
develop theoretical models including this new constraint. 



V. CONCLUSION 

Single file diffusion of electrostatically confined and in- 
teracting millimetric particles was studied. Their erratic 
motion conveyed through the mechanical shaking of the 
set-up was proved to be well described by the Langevin 
equation and thus to be Brownian. The characteristic 
parameters of the equation are well determined by fits of 
their radial mean square displacements. 

Long time behavior of the angular mean square dis- 
placement was explored. We exhibit a subdiffusive be- 
havior with a continuous slowing down that leads to a 
slower increase than the t^^"^ growth observed in colloidal 
systems. We suggest this might be a consequence of a non 
hard-core interaction in a periodic system. 



APPENDIX A: THE LANGEVIN APPROACH 
FOR A TRAPPED PARTICLE 

Experimental measurement of the time-dependent 
mean square displacement of a trapped particle requires 
to pay attention to the fact that the initial conditions 
of the measured trajectories are not equivalent, therefore 
the results are not exactly the same as the one usually 
presented in the theories. After a short summary of those 
usual results, we present what one should expect from an 
experimental determination such as ours. 

We consider a ID particle of mass m and coordinate 



x{t) trapped in a parabolic well of stiffness K and im- 
mersed in a thermal bath characterized by a friction 
force —ax and a rapidly fluctuating force F which obeys 
{Fit)) — 0, where the notation (} denotes an ensemble 
averaging. Newton's second law of motion gives : 



d X dx , , 

m—r- = —a— Kx + t it). 

dt-^ dt 



(Al) 



The major hypothesis of the Langevin equation is that 
F is fluctuating with a characteristic time Tc that is 
much lower than the time r^j = 7"^ that character- 
izes the variation of the speed of the particle, where 
7 = a/m is the damping coefficient [T^ . Then we can 
write {F(t)F{t')) = gS{t - t'). Note that 7 has the same 
microscopic origin as the strength F, that's to say the 
numerous interactions with the particles of the thermal 
bath, then we should expect 7 to depend also on g, there- 
fore on T. 

If one considers a given initial condition (xo,vo) and 
solves Eg. lAll one finds, in the case of a free particle (K 
—0), the following short and long time behaviors for the 
mean-square displacement Ax^(i) = {{x{t) — {x{t)))'^) : 



Ax^it) 



2Dt. 



(A2) 
(A3) 



where D 



is the diffusion constant. 



2771^7^ 7727 

For a trapped particle, the Langevin equation remains 
linear and is e asily s olved. In the case of weak friction 
(7 < 2 Wo = 2-^ K/m), one finds, for fixed (a;o,uo) initial 
conditions: 



Ax^{t) 



9 



0J~ 



-7* 



X (-^ cos{2Lut) - -Y sin(2tji) - ujI)\ , (A4) 



where = Wq — ^ . 

In particular, for long times {t th) 



Ax^(t) 



2m'^jujQ K 



(A5) 



However, in a majority of experiments, the averaging 
is obtained by taking different trajectories considered 
as different realizations of the random theoretical tra- 
jectory. Those trajectories are obtained by shifting the 
initial time, that's to say, if x{t) is the measured trajec- 
tory relatively to a given reference point, one considers 
the trajectories x{t-\-ta) — a;(io), where each value gives 
a new trajectory. 

For a free diffusion, this implies, since all points in 
space are equivalent, to take trajectories with the ini- 
tial condition = 0. However, initial speed cannot 
be controlled and summation over the trajectories will 
also implies an averaging over the distribution. If 
the experiment has started for a time larger than tj^. 
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we can consider that this distribution is the stationary 
one, thus the first moment of Vq is zero and the second 
is Note that this artefact cannot be avoided un- 

less one is able to start all experiments with exactly the 
same conditions. Denoting (()) this new averaging, we 
find ((x(t))) = and, noting that the vq and F distri- 
butions are independent, the following behavior for the 
mean square displacement : 



t-^o m ' 
- 2Dt. 



(A6) 
(A7) 



distribution. Again, if the experiment has started for a 
time larger than t^, we can consider that this distribu- 
tion is the stationary one, then that the first moment of 
xq is zero and the second is ^j^- 
We find that {{x{t) - xq)) = and 



Ax^{t) 
x[l- 

In particular : 



2((x^)) 



e ^ * ( COS! 



( cos(wt) + — sin(wt))] . (A8) 



The main consequence is that Ax'^{t) starts as t'^ for 
short times, and not as fEo. IA2p . Long time behavior 
is not modified. 

For a trapped diffusion, not all starting positions are 
equivalent and the averaging on the trajectories x{t) — xq 
necessarily implies averaging not only on the random 
force distribution but also on the initial condition {xq, vq) 



Ax' 



(A9) 



Short time behavior is thus independent from the well 
and similar to the free diffusion case. Note also that the 
infinite time limit ^^^^is twice the theoretical value in 
the case of identical initial condition trajectories. 
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